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We study the long time behavior of light particles, e.g. an electron swarm in which 

o> ' 

O 



Coulomb interactions are unimportant, subjected to an external field and elastic collisions 



with an inert neutral gas. The time evolution of the velocity and position distribution 

D. 

function is described by a linear Boltzmann equation (LBE). The small ratio of electron to 
neutral masses, e, makes the energy transfer between them very inefficient. We show that 
under suitable scalings the LBE reduces, in the limit e — > 0, to a formally exact equation for 
the speed (energy) and position distribution of the electrons which contains mixed spatial 
and speed derivatives. When the system is spatially homogeneous this equation reduces 
to and thus justifies, for e small enough, the commonly used "two-term" approximation. 



Introduction 

The motion of an electron under the influence of an external electric field and elastic 
"collisions" with the "background" is a problem of both historic and contemporary interest. 
It arises in gas discharges in the laboratory and in the atmosphere, when the degree of 
ionization is so low that electron-electron and electron-ion interactions are negligible, the 
so called swarm approximation [Carron 1992, Allis 1956, Franklin 1976, Margenau and 
Hartman 1948]. Similar situations arise in the study of electron transport in solids in a 
semi-classical regime [Lorentz 1916, Sommerfeld 1956, Ben Abdolllah and Degon 1996, Ben 
Abdolllah et al 1996, Golse and Poupaud 1992]. In such situations it is often reasonable 
to treat the scatterers as if they were in thermal equilibrium. For electrons moving in 
a (almost) neutral gas with constant density p this corresponds to the atoms having a 
Maxwellian velocity distribution F(V; M) at a specified temperature T, 



where M is the mass of the scatterers, e.g. the neutral atoms in a noble gas, k is the 
Boltzmann constant, and d is the space dimension. 

The time evolution of the electron distribution function is then given by a linear 
Boltzmann equation (LBE) of the form [Franklin 1976, Landau and Lifshits 1993], 



where e, m are the charge and mass of an electron, E is the external field and K gives the 
rate for scattering. The transition rate K will depend on the nature of the scatterers but, 
for the case of elastic collisions, to which we restrict ourselves here, K will satisfy detailed 
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(1) 



df(r,v,t) 



+ v— /- — ~/= [K(v,V)f(r,V,t)-K(v',v)f(r,v,t)}dV = JCf, (2) 
or m ov J 



dt 



balance with respect to a Maxwellian velocity distribution for the electrons F(v; m), at the 
same temperature as that of the neutrals, so that K(v,v')F(v';m) = K{y', v)F(v; m). 
This ensures that for E = 0, F(v; m) is a stationary velocity distribution which will be 
approached by /(v, t) as t — > oo; /(v, t) is obtained from /(r, v, t) by integrating over the 
spatial coordinate r. We assume for simplicity that the electrons are confined to a large 
periodic box: / is then a probability density with J J /(r, v; t)dvdv =1. It is also possible 
to consider other boundary conditions as well as electric and magnetic fields which are 
space and time dependent, but we shall not do so here. We shall also not consider here 
explicitly electrons in solids for which the scatterers are various kind of excitations [Ben 
Abdolllah and Degon 1996, Ben Abdolllah et al 1996, Golse and Poupaud 1992]. 

The central physical fact about the electron-neutral system is the great disparity in 



their masses: 



M 



e < 10 4 . Consequently the change in the speed of an electron in a 



typical collision, which is of order ^/e, will be very small, while the change in the velocity 
direction will be large, becoming independent of e when e — ► [Franklin 1976, Landau and 
Lifshits 1993]. This suggests approximating the integral operator in (2) by expanding the 
right side of (2) in powers of e and dropping "small terms" in e. There is a certain amount 
of ambiguity in carrying out such an expansion, arising from the uncertainty of how to 
treat the dependence of the unknown / on e. A reasonable choice gives 



df_ df_ _ eE df_ _ epd_ 
dt dr m dv v 2 dv 



v 4 a(v) fo + 



kTdfo 
mv dv 



+ vpa(v)(f -f). (3) 



In (3) v = |v|, /o( r 5 v, t) is the average of /(r, v, t) over angles and a{v) is the collision cross 
section, see [Margenau and Hartman 1948, Ginzburg and Gurevich I960]. The dependence 
of / on e is then to be determined from the solution of (3). 
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There is, however, very little control over this expansion. In fact eq.(3) is not necessar- 
ily positivity preserving. Nevertheless eq.(3) yields reasonable answers for the stationary 
velocity distribution of the electrons, /(v). The latter can be obtained explicitly in the 
so-called "two-term" approximation [Lorentz 1916, Margenau and Hartman 1948, Allis 
1956, Carron 1991], in which one keeps only the first two terms in a spherical harmonic 
expansion of /(v). This distribution was first found by Druyvesteyn [Druyvesteyn 1930, 
Druyvesteyn and Penning 1940], who considered the case T = 0, a = const, and was later 
generalized by Davydov [Davydov 1935] for all T and a. For an analysis of the error made 
by the two term approximation to (3) see [Rokhlenko and Lebowitz 1997, Rokhlenko - 
submitted] . 

A somewhat different approach to this problem was taken by Koura [Koura 1987]. 
Starting with a two term e-expanded kinetic equation (slightly different than what is 
obtained from (3)) he observed, that after a scaling of space, time and electric field and 
neglecting the time variation of the first harmonic, one is led to an equation in which e 
does not appear at all. Koura then argued that the actual physical quantities of interest in 
an electron swarm should have a similar simple scaling dependence on e, when e is small. 
This has the advantage of permitting more efficient computer simulations at larger than 
realistic value of e. Doing simulations for several values of e Koura found good agreement 
with results from the scaled two term approximation for e^lO -2 . 

In this note we shall use the same scaling as Koura but apply it directly to (2) without 
any other approximation. This is in the same spirit as the scalings used for electrons in 
solids [Ben Abdolllah and Degon 1996, Ben Abdolllah et al 1996, Golse and Poupaud 
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1992]. It is based on a formulation of space-time scalings now commonly used to obtain 
a "reduced hydrodynamic description" from a more detailed microscopic one, see [Spohn 
1995]. 

One-dimensional problem 

To simplify matters we shall first consider the one dimensional version of (2), corre- 
sponds to hard collisions of point particles with masses m and M = m/e, 

df(x,v,t) ^ df(x,v,t) _ eE df(x,v,t) = 
dt dx m dv 

(4) 

(l-e)V + 2ev m 



/oo 
-oo 



f(x,v,t)F(V;m/e)}dV, 



2V-(l-e)v 

T — t 

1 + e [ 1 + e e 

where vGl. Eq. (5) is to be solved subject to some initial condition f(x, v, 0). 

To obtain the behavior of / for long times we rescale our variables by setting y = ^fex, 
t = et (diffusive scaling [Spohn 1995]) and E = \ftE* (the field has to be small on this 
scale for the energy to remain bounded when the time is of order e -1 ). We now define the 
even and odd parts of the velocity distribution function 

1 

My, t) = 2^[/(y/-v/e, v, r/e) + f(y/y/e, -v, r/e)], 

(5) 

Vetp e (y, v, t) = ^=[f(y/^/e, v, t/e) - f(y/y/e, -v, t/e)], 

where v = |v] and we have put / — > 4j/ to preserve the normalization: in the scaled 
variable y the system is in a periodic box of size L independent of e. By changing the 
integration variable in (3) we then obtain two coupled equations for e and ip e , 

d<p e (y,v,r) dipt eE* dip e 

a + JT~ = 6 ) 

or ay m ov 



6-V 



m 



2nkT 



\v-V\ 



1 + e 
1 -e 



(l + e)v-2V 

y, ; ,t 

1 — e 



-My> v > T ) ) e 2ekT dv, 



dip e (y,v,r) dtp e eE* d<j> e 
e h v- 



m 



2ivkT 



\v - V\ 



di 



1 + e 
1 -e 



dy 



m dv 



(l + e)v-2V 

y, 1 ,t 

1 — e 



+ ip e (y, v, t) > e ™ckTdV. 



We assume now that the initial distribution is such that e and have well defined 
limits <j)(y, v, 0) and ip(y, v, 0) as e — > 0. Taking now formally the limit e->0on both sides 
of (7), keeping r, y and E* fixed we get the limiting equations, 



d<j)(y, v, r) ^ 9V» eE* dip ^ d 
dr dy m dv dv 



kT d(j> 
m dv 



d(j> eE* 

V-z 



-2pvip(y,v,T), 



(8a) 
(86) 



dy m dv 

valid for v > 0. Solving (8b) for ip and substituting into (8a) we get the reduced equation 
for the even part of the distribution (f)(y,v,r), in terms of the scaled space, time and 
electric field 



dcj) 1 d (eE* d(p _ v d<f>\ __ d 
dr p dy \ m dv 2 dy J dv 



with the boundary condition at v = 0, 



eE 



1 kT 
— + —2pv 
m J Zpv m 



tv^. 



(9) 



dc/>(y,0,T) _ eE* Um ld(j>(y,v,T) 



dy 



m v^+o v dv 



(10) 



and the normalization 



I dy I (f)(y,v,r)dv = 1/2. 
Jo Jo 



The condition (10) at v = follows from the (assumed) continuity of the distribution / at 
v = 0. 
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Let us introduce the rescaled electron density n(y, r), drift u(y,r), and mean speed 
w(y,r): 

POO 

n(y,r) = 2 c/)(y,v,T)dv, 

J o (11) 

pOO poo 

n(y,r)u(y,T) = 2 vi(;(y,v,T)dv, n(y, r)w(y, r) = 2 / v(f)(y,v,T)dv. 
Jo Jo 

Integrating Eq. (9) over velocities, we obtain the continuity equations 

dn(y,r) d{nu) = Q 



Or 



dy 



(12) 



To obtain u we substitute (8b) into (11) and integrate over v, 



e0(y,O,r) w(y,r) dn 1 dw(y, r) 

U = ; zti 



(13) 



mpn(y,T) 2pn(y 7 r) dy 2p dy 
One can identify in (13) the factors in front of E*, ^ as respectively the mobility, 
diffusion coefficient, and a parameter related to thermodiffusivity [Golant et al. 1980]. 

The right side of eq.(9) coresponds to a drift-diffusion in "speed space" with a diffusion 
coefficient given by the term in the square brackets and with a drift —2pv 2 . The spatially 
hompgeneous stationary solution of (9) has the form 



4>{v) = C exp 



Jo 



s 3 ds 



MLo2 i ( eE* V 
m° ^ \ 2mp J 



(14) 



similar to the Davydov distribution for d = 3, c.f. eq. (21). Note that for E* = 0, i.e. no 
external field, 4>(v) = F(v; m) as it should. 



Three dimensions 



In d = 3 the right side of the LBE takes the form 



Kf = pj J |v - VKv - V) [f(V, t)F(V; M) - /(v, t)F(V; M)] JVdu 



The angular integration is over the scattering solid angle Cj of the electron in the center 
of mass coordinate system, cr(v) is the electron- neutral differential collision cross section, 
and we have 



/ = v+ — - (V-v)], 



2e 



V = V Cv[u ■ (V-v)]. 

1 + e 



Carrying out a similar analysis as for d = 1 we separate / into a spherically symmetric 
part e (y, u, r) and a remainder ^ e (y, v, r), i.e. we set 



/(y/v^, v, t/e) = e 3/2 [(/> e (y, u, r) + v^My, v, r)], 



e -3/2 /• 

^e(y,v,T) = -^ r - / f(y/y/e, v,t/e)(ffi. 



(15) 



(16) 



In (16) the integration is over the unit sphere specifying the orientation of the vector ^. 
Taking now formally the limit e^Owe obtain, in terms of the rescaled variables, a set of 
equations for <fi and ip entirely analogous to (8), (9), 



<9r 3/90" y 3mp 



1 90 1 9 /'w 2 ' 



i__a 

V 2 <91> 



v dcj) 



a dv v 2 dv \ a 
kT d<f) 



+ pv a(v) 4> + 



mv dv 



(17) 



m J 3ap dv 

iP = —L- (-E* • V v - v • V y </>) . (18) 
pa(v)v \m J 

We assume the collisions to be spherically symmetric, so a(v) = 47ra(v) is the total cross 



section. 



The spatial density and current are given by 



n(y,r)=47r/ v 0(y, v, r)dv 7 n(y, r)u(y, r). 

'0 



(19) 



They satisfy equation of continuity 



dn(y,r) 
dr 



+ V y [mi] = 0. 



(20) 



The electron drift u has the form 



u = uE* - D-V v n V 

n 3p 



y 



( 




i d r v 2 ' 



D = — 



H = 



3pm \v 



2 dv a(v) 



where (g(v)) = ^ / °° gv 2 4>{y, v , r)dv is the average over the velocity distribution. In the 
spatially homogeneous case eq. (17) coincides with the usual equation obtained in the two 
term approximation, whose stationary solution is 



which was given by Davydov in [Davydov 1935]. For E* = 4>(v) in (21) is just F(y;m) 
while for T = 0, a = const it coincides with the Druyvesteyn distribution. 

Inequalities for moments at arbitrary e and T=0 

In previous works [Rokhlenko and Lebowitz 1997, Rokhlenko - submitted] we estab- 
lished bounds on the first few moments of the stationary distribution function using the 
approximate equation (3) with a velocity independent total cross section a . Here we shall 
do the same for the LBE (2) when T = 0. Introducing the dimensionless variables, 



and noting that for T = 0, F(V;M) = 6(V), the stationary, spatially homogeneous, / 
satisfies, in d = 3, the equation, 




(21) 
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(1 + e) 2 



J f(s')d[(l - e)(s') 2 + 2e(s ■ s') - (1 + e)s 2 ]ds' - sf(s). 



(22) 



2tt 

Following [Lorentz 1916, Margenau and Hartman 1948, Rokhlenko 1991] we expand 
the distribution function in Legendre series 



f(8,x) = J2f l (8)Pi(x), 



1=0 

substitute it into (22) and obtain a set of coupled equations for the //. Let us introduce 
the moments 

poo 

Mi(k) = / s^^ds. (23) 



o 



We can now use the same technique as in [Rokhlenko and Lebowitz 1997] to obtain bounds 
on the mean energy W(e) and drift u{e) which are defined in terms of M (k). These 
inequalities are based on the fact that log M (k) is a convex function of k. This yield in 
the present case 

Vb ^ 2pa^ Tjr/ , ^ 1 6 1 / 4 ^ 2e 3 / 4 , s ^ 1 

-^^f w(e ^7S' T ~ w^V^ M(£) -^' (24) 



where 



= p (4,e)pi(2,e) p (7, e)pi(5, e) , , 

a 3e ' 12e[l+po(3,e)/p2(3,e)] 1 J 



and 



2 /■! 



Ze J\l-e/l+e\ 



(1 + e)x 2 + e - 1 



cfe - 1. (26) 



2ex 

When e — > oo the lower and upper bounds for u and W merge giving in the limit 



eE mu 2 , , 

u = J— , W=— — . 27 
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When e — > 0, a — > 1, 6 — > 1/2 and the inequalities (24) are then satisfied by the 
Druyvesteyn distribution. We believe that with greater effort it should be possible to 
obtain upper and lower bounds which when e — > both converge to the values obtained 
from the Davydov distribution. 
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